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ABSTRACT 

Current theories of spiral and bar structure predict a variety of pattern speed behaviors, calling 
for detailed, direct measurement of the radial variation of pattern speeds. Our recently developed 
Radial Tremaine- Weinberg (TWR) method allows this goal to be achieved for the first time. Here we 
present TWR spiral pattern speed estimates for M101, IC 342, NGC 3938 and NGC 3344 in order to 
investigate whether spiral structure is steady or winding, whether spirals are described by multiple 
pattern speeds, and the relation between bar and spiral speeds. Where possible, we interpret our 
pattern speeds estimates according to the resonance radii associated with each (established with the 
disk angular rotation), and compare these to previous determinations. By analyzing the high-quality 
HI and CO data cubes available for these galaxies, we show that it is possible to determine directly 
multiple pattern speeds within these systems, and hence identify the characteristic signatures of the 
processes that drive the spiral structure. Even this small sample of galaxies reveals a surprisingly 
complex taxonomy, with the first direct evidence for the presence of resonant coupling of multiple 
patterns found in some systems, and the measurement of a simple single pattern speed in others. 
Overall, this study demonstrates that we are now in a position to uncover more of the apparently 
complex physics that lies behind spiral structure. 

Subject headings: galaxies: individual (M101, IC 342, NGC 3938, NGC3344) - galaxies: spiral - 
galaxies: kinematics and dynamics - galaxies: structure - methods: numerical 



1. INTRODUCTION 

The angular rate at which bar and spiral structures 
rotate, the pattern speed, is a parameter appreciated for 
its capacity to reveal information about the appearance, 
properties and evolution of galaxy disks. Since the pat- 
tern speed determines the rate at which gas and stars 
encounter these structures, it provides a description of 
their influence on star formation, the processing of the 
interstellar medium and the dis t ribution of metals (e.g. 
iRancfl [19981 : iKnanen et afl[i"996l iHenrv et al.l[2003h . for 
example, as well as the redistribution of mass in the disk. 
Bar torques dep end on the pattern speed, and therefore 
gas inflow rates (jQuillen et al.lll995l ) do, as well. In early- 
type barred galaxies, pattern speeds have also been used 
to investigate the c entral Dark Matter (DM) content in 
these systems (e.g. Debattista fc Sellwoodll2000t iCorsinil 
[20081: but see lDubinski et al.ll2009l ). 

The pattern speed moreover supplies a fundamental 
characterization with which to investigate the origin and 
evolution of large-scale structures, one of the prime un- 
resolved issues in the dynamics and evolution of galaxy 
disks. With knowledge of bar and spiral or the radial 
variation of spiral pattern speeds in a given galaxy it is 
possible to determine whether spiral structure is steady 
or winding, the domain and number of patterns that 
can be sustained in a disk, and the relation between bar 
and spiral pattern speeds. Direct pattern speed measure- 
ments, however, have been elusive. 



Although grand-design spirals have been observa- 
tional ly linked to the pre s ence of bars or compan- 
ions (jKormendv &: Norman! Il979t ). it remains unclear 
whether spirals are long-l ived density waves persisting 
over many revolutions (i.e. iLin &: Shu 1964; Bcrti rTet al.l 
1989), or rapidly-evolving transient features, as 
found in several simulations ( Scllwood & Carlber gj ll984h 
iThomasson et a.lJfl99(l : ISellwood fc Kahnlll991h . 

It is also uncertain how waves can exist over large frac- 
tions of a disk radius, as observed. Since the resonances 
(inner and outer Lindblad; ILR and OLR) between which 
waves are expected to propagate in general do not span a 
large range in radius, waves of different speeds and struc- 
ture m ay occupy distinct radial zones (|Masset fc Taggerl 
Il997bf ). In barred galaxies, the apparent alignment be- 
tween bar and spiral (as in M83, where the two-armed 
spiral emanates from the bar ends), initially taken as ev- 
idence f or a common pattern spe ed, does not occur in 
general (jSellwood fc Sparkelll988f h Furthermore, where 
the bar ends near its Corotation Re sonance (CR) (e.g. 
as found in early-types: ICors ini 2008), a spiral generated 
with the same speed would lie mostly outside the CR 
(i.e. with dust lanes along the convex side of the arms), 
which is also uncommon. 

In the theory of "mod e coupling" (Sv gnet et el! .1119881 
iMasset fc Taggerl Il997bl ). multiple patterns in different 
radial zones are linked, such that the resulting wave 
structure can extend over a larger radial range than 
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is possible for a single pattern. In this scenario, the 
CR of an inner pattern overlaps with the ILR of an 
outer pattern, and at this overlap energy and angu- 
lar momentum are efficiently transferred outward in the 
disk. This was fir st demonstrated between simulated 
ba rs and spirals bv iMasset fc Taggerl (|1997b) and later 
bv lRautiainen fc Said (jl9991 who find, in addition, that 
other resonance overlaps, such as the coincidence of the 
spiral's inner 4:1 resonance with the bar CR, are also 
effective. Given that the radial extent of resonant or- 
bits can be quite large, these overlappings can be spa- 
tially broad, and frequency diagrams can therein be suit- 
ably revealing. However, not all such resonance over- 
laps are accompanied by signs of true mode-coupling (i.e. 
boosted beat mod es detectable in simulation power spec- 
tra at the overlap; [Masset fc Taggerl (|1997bf )). 

Spiral-spiral mode coupling may also occur, typified, 
perhaps, by the transition commonly observed between 
a strong two- armed spiral and more complex structure 
at som e radius. In their simulations, iRautiainen fc Said 
( 1999) find evidence for spiral structure in the absence of 
a bar, spiral-spiral mode coupling, and multiple pattern 
speeds without mode coupling. 

To observationally investigate these possibilities, 
we have undertaken a project to measure spi- 
ral pattern speeds and their radial variation using 
the Radial Tremaine- We in berg method developed by 
iMerrifield. Rand fc MeTdtl (120061) . su bsequently tested 
on simulations (jMeidt et alJ l2008af ): hereafter Pa- 
per I), and recently applie d to a real spiral in the 
grand-design galaxy M51 (jMeidt et all l2008bl here- 
after Paper II). Like its traditional counterpart-the 
TW m ethod for measuring a sing le, constant pattern 
speed (|Tremaine fc Weinberg] Il984h -the TWR calcula- 
tion yields measurements of pattern speeds using ob- 
servationally accessible quantities based on a require- 
ment of continuity. The method is direct, and so 
overcomes several of the obstacles in reliably estimat- 
ing pattern speeds with other, indire ct methods (e.g. 
reson a nce identification or modeling; [E lmcgrcc n et all 
19891 lElmegreen et all 119961 IRautiainen fc Said 120051 . 
Garc ia-Burill o et al.lll993 



Moreover, the TWR method provides a resolution for 
those applications of the TW method where pattern 
speed estimates sh ow a clear, systematic departure f rom 
a single value (e.g. IZimmer. Rand fc McGrawll2004l and 
IMerrifield. Rand fc Meidtl I2006D . Particularly in disks 
with multiple or extended structures, this behavior im- 
plies that the pattern speed is not constant either be- 
cause it varies temporally or spatially. By allowing the 
pattern speed to vary in the radial direction, the TWR 
method explicitly allows for the possibility that a galaxy 
may contain a number of distinct features at different 
radii, such as bars and spiral arms, each with their own 
pattern speeds. It also makes possible the detection of 
spiral winding and hence the estimation of the lifetime 
of a galaxy's current pattern. 

Our aim with this paper is to expand to a larger sam- 
ple of galaxies the TWR analysis recently applied to 
the grand-design spiral galaxy M51 (Paper II). Here, as 
there, we apply the method according to the prescrip- 
tion outlined in Paper I. The calculation employs reg- 
ularization, which smooths otherwise intrinsically noisy 
solutions through the use of a prior models of the radial 



dependence of the pattern speed, and affords straightfor- 
ward tests for bar-spiral and spiral-spiral relations, and 
spiral winding. 

Also as in Paper II, for use as a kinematic tracer we 
consider observations of the ISM, which have become the 
standard choice of spiral tracer for meeting the conti- 
nuity requirement of the method. The application of 
the TW and TWR met hods to CO and HI observations 
of sp i ral galaxies (e.g. Wcstpfa hlll998t iRand fc Wallinl 
12004 IMerrifield. Rand fc Meidtl I2006D avoids significant 
problems with the stellar component in these systems, 
namely the faintness of the old stellar disk, and the ef- 
fects of star formation and obscuration by dust in spiral 
arms by which the application of the continuity equation 
is invalidated. 

Furthermore, where the ISM is dominated by either 
the molecular (traced by CO) or atomic gas phases, 
conversion among phases can be assumed to occur at 
low levels on orbital timescales such that, together with 
the low true efficiency of star formation in spirals, the 
dominant component arguably obeys continuity (e.g . 
IZimmer. Rand fc McGrawl l200l IRand fc Wam^Bool . 
In addition, as an improvement on previous ISM-based 
TW spiral studies, here we analyze both molecular gas 
observations from the BIMA SONG (|Helfer et all 12001 
maps include zero-spacing flux information) and archival 
21-cm emission data tracing the atomic hydrogen phase. 
(The flux information at the largest scales in all of the 
HI cubes considered here is comparable to that achieved 
with single-dish observations; see references in § §2.2112.5]) . 
Our consideration of both CO and the more extended 
HI here serves two main purposes. In galaxies that are 
not molecule-dominated over the extent of the detectable 
CO emission, the HI supplements H2 to establish a total 
particle, continuity-obeying tracer. The radial range of 
detectable pattern speeds is also increased in this way. 

Part of our treatment, in this case, entails an investi- 
gation (where possible) into the sensitivity of TWR so- 
lutions to the CO-to-H2 conversion factor X adopted in 
combining the data, which h as been suggested t o vary lin- 
earl y with metallicity (e.g. iBoselli et all l2002h . Follow- 
ing IZimmer. Rand fc McGrawl (|2004f ) and Paper II, we 
consider the effect of variation in X with radius (but not 
the possibility, for example, of arm-interarm variations). 
We also take into consideration distortions or warps ob- 
servable in the outer HI disks in our sample (described in 
§ §2.1112.5]) . which violate of one of the main TW assump- 
tions (namely, that the disk is flat: lTremaine fc Weinberg! 
fl98l . 

Based on the a priori models described in §3.11 the 
best-fit pattern speed solutions calculated for each galaxy 
are presented in §§ 13.2113.51 There, the resonances asso- 
ciated with each, which we identify through comparison 
with angular rotation curves, serve as a main informant 
of our results; our findings are motivated by, and com- 
pared with, the large body of work linking resonances to 
the dynamics and morphology of disks. We also discuss 
the limitations and sensitivities of the TWR calculation 
(described in Paper I) as applied to each galaxy and sum- 
marize our results in §4. 



2. THE SAMPLE 
2.1. Selection and Overview 
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Fig. 1. — Clockwise from top left: Zeroth moment maps for M101, IC 342, NGC 3938 (made from the combined CO and HI cubes) and 
NGC 3344 (made from the HI cube) showing the logarithm of the column density N in units of cm -2 . The horizontal bar near the bottom 
indicates the physical scale. 



In the following subsections we describe the four spi- 
ral galaxies analyzed in this paper, M101, IC 342, NGC 
3938, and NGC 3344. This small sample is not meant 
to be representative of variations in spiral structure with 
galaxy classification type, or to embody all possible the- 
ories of spiral structure. These nearby galaxies were se- 
lected simply for their clear spiral structure and mod- 
erate inclination, as well as for the availability of ob- 
servations (specifically, both CO and HI in three of the 
four cases) of sufficient sensitivity and resolution to allow 
structure to be resolved and the radial variation of Q p to 
be investigated; see Paper I. 

With this sample we explore three scenarios for ap- 
plying the TWR method to observations of the ISM. We 
first consider both the molecular and atomic components 
of the ISM in M101, which is molecule-dominated over 
the detectable extent of the CO emission, and investi- 
gate whether there exists a relation between the speeds 
of the multiple structures apparent throughout the disk. 
We then apply the method to two galaxies (IC 342 and 
NGC 3938) where the ISMs are dominated by neither 
the atomic or molecular gas phases in the central few 



kpc; in this case, considering the total column density is 
necessary for meeting the continuity requirement of the 
method. Lastly, for NGC 3344, where the ISM is domi- 
nated by the atomic phase over the majority of the disk 
area (see §2.5|) . we analyze HI data alone, and with our 
TWR solution predict a pattern speed for the small bar. 

In three of the four cases, a warp-like distortion is ev- 
ident in the HI distribution and kinematics, often with 
little further evidence of a strong spiral pattern. A warp 
in the outer disk, characterized by position and/or in- 
clination angles that differ from those values in the disk 
interior, violates the TW assu mption that the disk is flat 
(jTremaine fc Weinberg! [l984f ). In implementing the reg- 
ularized TWR method we therefore adopt the procedure 
advocated in Paper I (see the beginning of § I3.ip , which 
accommodates for the presence of information about the 
warp in the TWR quadrature by excluding it from reg- 
ularized solutions. 

Zeroth and first moment maps for each galaxy are 
shown in Figures Q] and H For M101, IC 342 and NGC 
3938, these maps were generated from the combined CO 
and HI cubes, which were regridded to a uniform channel 



4 




Fig. 2. — Clockwise from top left: First moment maps for M101, NGC 3938, IC 342 (made from the combined CO and HI cubes) and 
NGC 3344 (made from the HI cube) in units of km s —1 . 



width and co-added in units of particles cm~ 2 . Prior to 
combination, the CO cubes were smoothed to the resolu- 
tion of the HI and CO intensities converted to molecular 
column densities assuming a constant CO-to-i?2 conver- 
sion factor X=2. 0xl0 20 cm- 2 [K km s -1 ]" 1 , the Galactic 
mean value (e.g. iHunter et al.1 [19971. (Subsequent vari- 
ations in the X-factor are considered when applicable.) 

Tilted ring fits to the velocity fiel d of the four galaxies 
were performed using the GIPSY (jvan der Hulst et all 
1992) task ROTCUR in order to derive the kinematic 
parameters used in the TWR calculation. For these fits, 
the systemic velocity is initially fixed to the value given 
in the literature, while the kinematic center, position an- 
gle (PA), and inclination are allowed to vary. (The value 
and the errors on V sys were subsequently determined by 
fitting with all other best-fit parameters held fixed.) The 
values we find are in good agreement with those nomi- 
nally adopted for each galaxy and are listed in Table [TJ 

For the purposes of identifying resonances (corotation, 
inner and outer Lindblad), we also derived the circu- 
lar velocity in each ring with the best-fit kinematic pa- 
rameters held fixed. In each galaxy these velocities are 
modeled in the unwarped portion of the disk with least- 



squares fits to the standard three- parameter approxima- 
tion (i.e. iFaber fc Gallaghenll97ll 



V ro t(r) 



Vrnax y / ^max ) 



(i) 



(l/3 + 2/3(r/r Wc )«) 3/2n 

as shown in Figure [3] The errors bars on each mea- 
sured velocity there reflect the average deviation from 
this value on either the approaching or receding side, 
which together we find better represent the uncertainty 
in the rotation curve than do the formal errors returned 
by ROTCUR. With the fitted velocities, we then gener- 
ate a set of smooth curves for fi, fi±«/2, and J7±k/4; 
these curves are intended to reduce the impact of non- 
axisymmetric motions (e.g. spiral streaming) on our res- 
onance identifications and were invoked without regard 
to specific mass component characterization. 

2.2. M101 

This SABcd galaxy (D=7.4 Mpc; Uurcevic fc Butcher! 



HH) is an excellent candidate for analysis with the 
TWR method. The molecular gas filling the central 3' 
hole in the HI emission f e atures a bar discovered by 
iKennev. Scoville fc Wilsor] (|1991f h The bar is a clear 
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Fig. 3.— Clockwise from top left: Rotational velocities for M101, IC 342, NGC 3938, and NGC 3344. Values measured with tilted ring 
fits to the velocity fields shown in Figure [2] with all optimal parameters listed in Table [T] held fixed are shown as filled circles (crosses) in 
the unwarped (warped) region of the disk (as identified in §§ 13.21331 . The points are shown spaced at the resolution of each map with 
error bars (shown in dark gray) representing the average variation in the derived velocity from side to side (approaching and receding). 
The solid line in the plots for M101, IC 342 and NGC 3938 is the least-squares fit of the velocity model (Equation [TJ to the velocities 
measured in the unwarped portion of the disk (filled circles), with best-fit parameters: Vmax=189 km s — , rma X =15.89 kpc, n=0.43 for 
M101; V max =167 km s" 1 , r max =5.62 kpc, n=0.76 for IC 342; y ma:c =154 km s" 1 , r max =7.3 kpc, n=0.76 for NGC 3938; V max =l68 km 
s , r max =5.8 kpc, n=0.29 for NGC 3344. (For NGC 3344, to achieve an accurate fit these parameters had to be constrained a priori to a 
range that best reprodu ces the flatness for r>3 kpc). All solid lines are used to derive the smoothed angular rotation and resonance curves 
presented in §§ 13.213.51 



TABLE 1 

Parameters used in the TWR calculation. 



Parameter 


M101 


IC 342 


NGC 3938 


NGC 3344 


Dynamical Center RA (a) (J2000) 


14 ft 3 m 13M3 


3 h 46 m 48 a .4 


ll A 52 m 49*.8 


10 ft 43 m 31 a .5 


Dynamical Center DEC (5) (J2000) 


54°20'56" 


68°5'47".8 


44°7'11".7 


24°55T8".3 


Distance (Mpc) 


7.4 


2.0 


11.3 


6.9 


Systemic Velocity (Vsys, km s _1 ) 


244±8 


30±2 


809±3 


586±3 


Position Angle(°) 


42±3 


42±3 


21±2 


155±2 


Inclination (°) 


21±6 


31±5 


M±3 


25±4 



Note. — Entries for the dynamical center, inclination and position angle for each galaxy were derived with a tilted ring analys is of the 
first moment of the data cube using the GIPSY task ROTCUR. Optimal values for the distance originate with references cited in 12 . 212 ."51 



candidate-driver of the spiral structure that appears to 
emanate from the bar region. At the outer radii, tidal 
interaction wit h companion galaxies (e.g. NGC 5474 
andNGC 5477: iHuchtmeier fc Witzelfl97a IWaiTer et all 
1997) is thought to be responsible for the lop-sidedness 
and distortion in th e HI distribution, a nd may also be the 
source of the spiral ([Waller et al.lll997l ) . The high resolu- 
tion (~7") of the total H2+HI maps constructed from the 
BIMA SONG CO data a nd THINGS HI dat a (with RO- 
BUST weighting scheme: [Walter et al.ll2008h allows us to 
apply the TWR method with exceptional leverage on the 



radial dependence of the pattern speed throughout the 
disk. We aim to derive the bar and spiral pattern speeds, 
and identify multiple spiral modes and spiral winding, if 
present. 

As considered in the upcoming discussion f§ 13. 2[) . we 
note that the rotational velocities between r~7 and 14 
kpc and at radii r>19 kpc in the top right of Figure[3]are 
not well fit here. Also, while the rise in the latter zone 
appears on both the approaching and receding sides, a 
^75 km s" 1 asymmetry exists the rein. 

In addition, as demonstrated by IWong fc Blitzl ([2002) 
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assuming X=1.8xl0 20 cm _2 [K km s -1 ] -1 (slightly lower 
than our adopted value), and using the BIMA SONG 
data we analyze here, the ISM is molecule-dominated 
over the extent of the CO emission. However, the disk 
is known to sustain a metallicity gradient, which, for a 
linea r dependence of X on metallicity (e.g. iBoselli et alJ 
2002), could imply variation in X with radius. So, while 
a constant offset should have no implications for TWR 
solutions, we also consider the effect of variation in X 
outward from the center. 

2.3. ICS42 

For this Seyfe rt SABcd galaxy (D=2.0 Mpc; 
ICrosthwaite et all [2000), we base our analysis on the 
BIMA SONG CO data cu be and the 38" VLA HI o bser- 
vations first published bv ICrosthwaite et alJ (|2000f ). As 
in M101, the CO fills the central hole in HI distribution, 
but here the atomic gas makes a significant contribution 
to the total gas density before the edge of the CO emis- 
sion. A thorough discussion of the HI-CO overlap and 
the features (a bar, a two-armed spiral, and a four-armed 
spiral) in th e two gaseous componen ts can be found in 
the study of ICrosthwaite et alJ (|2001h . 

2.4. NGC 3938 

NGC 3938 (D=11.3 Mpc: Uimenez- Vicente et aflllQQfll) 
is a nearly face-on, late type (SAc) galaxy exhibiting a 
central two-armed spiral that branches into multi-armed 
structure in the optical. The two strong spiral arms are 
evident in the molecular gas, which reaches a surface den- 
sity comparable to that of the HI inside the edge of the 
CO emission. But where the stellar disk exhibits floc- 
culen t, but clear, spiral stru cture (with as many as six 
arms;lElmcgrc en" et al.lll992T ). the HI disk is characterized 
by less well-organized, irregular structure. The combina- 
tion of archival WSRT HI data with BIMA SONG CO 
observations establishes a total H2+HI kinematic tracer 
at roughly 20" resolution that e xtends to just beyond the 
optic al extent of this galaxy (|van der Kruit &: Shostakl 
[l982TV 

2.5. NGC 3344 

The application of the model-independent TWR 
method to thi s Hi-dominated, isolated SABb c galaxy 
(D=7.4 Mpc; IVerdes-Montenegro et~aTl I2000D should, 
in principle, clearly establish the relation of the ring- 
like morpholo gical features at r=l kp c and r=7 kpc, 
identified by IVerdes-Montenegro et alJ (|2000l ) in opti- 
cal images, to the large-scale patterns present in the 
disk (and their resonances). The outer ring, in par- 
ticular, although thought unlikely to be related to the 
small bar (not observed in the gas) with as^0.7 kpc 
(|Verdes-Montenegro et al. 2000), has not been otherwise 
conclusively related to the spiral. In search, as well, of 
spiral winding and multiple spiral modes, here we an- 
alyze 20" resolution archival WSRT HI data where the 
dominant two-armed spiral and the outer ring are clear, 
in addition to the outer lopsided region e xhibiting a twist 
in the isovels (Ivan der Kruit fc Shostaklfl98l . The zone 
of the bar and the inner ring, which falls within the cen- 
tral 26" where there is little 21-cm emission, corresponds 
to less than two resol ution elements. S o although the 
ring is resolved in CO (Re gan et al J 12002). and the peak 



H2 column density (|Helfer et al.ll2003[ ) there exceeds the 
HI, we do not consider the contribution of the molecular 
gas here. 

3. RESULTS 
3.1. Applying the TWR Method 

We apply the regularized TWR method as in Paper I. 
For each galaxy we consider several smoothed, testable 
models for f2 p (r). These models vary as polynomials (or- 
der n<2) designated into at most three distinct radial 
zones. Where a priori evidence suggests that there is 
little information from a strong pattern beyond a cer- 
tain radius, or that the TW assumptions are otherwise 
violated by the presence of a warp, for example, our mod- 
els also include the parameterization of a cut radius, r c , 
beyond which all bins are calculated without regulariza- 
tion (i.e. the functional form is unconstrained). These 
models, with r c marking the end of the dominant struc- 
ture, have been demonstrated to sufficiently separate the 
compromised zones in the disk from regions where infor- 
mation about patterns can be reliably extracted in the 
TWR calculation (Papers I and II) . All transitions r t be- 
tween distinct zones and all cut radii, where present, are 
treated as free parameters. 

For each model, the two numerical solutions for 
O p (r) from either side of the galaxy (y>0 and y<0; 
iMerrifield. R and fc Meidti [2006) ' are averaged to con- 
struct a single, global model solution. Each model solu- 
tion is then judged based on a simple reduced x 2 statis- 
tic, with the best-fit solution corresponding to the % 2 - 
minimum in the full parameter space. 

As in Paper II, the random, measurement errors used 
in the regularized calculation, and with which we judge 
the best-fit solution, reflect uncertainty arising with the 
chosen flux cut-off in the first moment maps. The sys- 
tematic errors on each measurement represent uncer- 
tainty in the PA, which is the d ominant source o f er- 
ror in TW and TWR estimates (pebattistal 120031 : Pa- 
per I), roughly 20% for 5pa=3°. But, here we report 
these as a dispersion on each measured value, rather than 
present individual solutions for each PA; this is possible 
here since, unlike in M51 (Paper II), we find no mean- 
ingful evidence that the form of the model associated 
with the best-fit solution for any of the galaxies in the 
current sample varies from PA to PA. Also, unless oth- 
erwise noted, errors due to uncertainty in, e.g., the incli- 
nation angle are generally smaller and are not reported; 
these prove to be of little consequence to the accurate 
placement of radial bins defined in the quadrature (as 
suggested in Paper I). The additional change introduced 
in the measurements Q p through a change in sin i is fur- 
thermore shared by £1 and k, and so our resonance iden- 
tifications, in particular, should not be effected by error 
in the inclination to first order. A thorough account of 
our methodology can be found in Paper I (and references 
therein) . 

3.2. M101 

We apply the TWR calculation with radial bin width 
Ar=7"=0.27 kpc (D=7.4 Mpc), the resolution of the 
combined cube. Together with the position of the outer- 
most slice on each side, | j/ 1 =30.4 cos i kpc, this establishes 
the extent of integration along each slice, equipping solu- 
tions with 113 bins in total. The best-fit solution given a 
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Fig. 4. — Best-fit regularized solution for M101 with 
r c =21.9±0.43 kpc for PA=42°±3°. For this solution, bins exterior 
to rc (not shown) have been calculated without regularization. 
Dashed red lines represent the dispersion in solutions derived with 
a three-pattern speed model at PA=39° and 45° . Horizontal error 
bars represent the dispersion in rt y i, rt,2 and r c from PA to PA. 
The innermost speed corresponds to n p i=47±10 km s _1 kpc _1 
out to rt, i=6. 7±0. 25 kpc, followed by Q P: 2~18±1 km s _1 kpc _1 
out to ^^=13.8^10.58 kpc and H Pi 3=5±3 km s — ^pc" 1 out to 
tq- Curves for Q, Q±k/2 and £2±k/4 (see text) are shown in 
black, cyan and blue. 
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Fig. 5. — Fourier power spectrum of the zeroth moment map 
shown in the top left of Figure [l] Modes up to m=4 are plotted as 
a function of radius with lines for m=l in black dash-dot-dot-dot, 
m=2 in red solid, m=3 in green dash-dot, and m=4 in blue dash. 

PA uncertainty ±3° is plotted in Figure [4] There shown, 
also, are the rotation and resonance curves derived as 
described in ij2.ll It should be noted that a ~75 km 
s _1 asymmetry in the rotational velocities from the ap- 
proa ching and receding sides e xists at radii r>20 kpc 
(e.g. IKamphuisl[l99a IJodl2002l) . 

The outermost portion of the disk is effectively re- 
moved from the solution with the parameterization of 
a cut radius r c =21.9±0.43 kpc. This radius, identified 
at the minimum of the x 2 , is comparable to the loca- 
tion where the disk becomes visibly distorted. We find 
that the outer distortion/lop-sidedness clear in the sur- 
face density is well characterized by the predominance 
of an m=l asymmetry beyond r~20 kpc in the Fourier 
decomposition (Figured]) and also matched by a warp in 
the outer velocity field; in addition to the rotation curve 
asymmetry, with our ROTCUR analysis we find that the 
PA and inclination of fitted rings beyond this radius vary 
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Fig. 6. — Total H2+HI surface density in M101 highlighting the 
structure inside ~18 kpc. The transitions rt 1 and rt 2 marking the 
extents of the first and second pattern speeds in our solution are 
shown as dashed white circles. The corotation radius rcii=3.5 of 
the innermost speed is shown as a solid white circle. The horizontal 
bar near the bottom left indicates the physical scale. 

significantly from the nominal values established in rings 
interior. (We note that this strong variation in the PA is 
found to start at a much large r radius than identified by 
iRownd. Dickey fc~ Hclou 1994.) Inclination variation, in 
particular, is a violation of the TW/TWR assumptions 
and so we argue that, by excluding the bins covering the 
outer disk from models for Q p (r), the remaining, inner 
regularized bins are better equipped to reproduce the 
true pattern speed. (Note that with this cut radius a 
possible m—1 mode describing the outer, lop-sided por- 
tion of the disk is ignored). 

In this case, the most conspicuous aspect of the so- 
lution is the pair of transitions at r tj i=6.7±0.25 kpc 
and r t 2 = 13.8±0.58 kpc (marked in Figure [H]) between 
three distinct pattern speeds, f2 Pi i=47±10 km s _1 kpc _1 
f2 P) 2=18±l km s _1 kpc _1 and Pi 3=5±3 km s _1 kpc _1 
out to r c . As discussed below, the innermost, con- 
stant speed is nearly identical to the value f2 p =48±8 
km s~ kpc - measured from the CO maps, alone (see 
§ 13.2. ip . while the second pattern speed is consistent 
with the value generally upheld in the literature (O p ~19 
km s _1 kpc _1 for r C R~12 kpc, e.g. IWaller et all 119971 ; 
but see the end of this section) . 

The quality of our measurement is especially clear in 
the precision with which the transitions between the 
three pattern speeds are determined. The transition r tj i, 
for example, shows less than 5% variation from PA to 
PA. This transition coincides with a marked decrease in 
m=2 power (Figure [5]) and occurs well past the edge of 
the CO emission. We can therefore recognize that the in- 
ner speed describes the molecular bar and the two-armed 
spiral manifest in the CO an d weakly traceable in HI, 
which Elmegrc en et al.l (|1992f ) identify in the B-band out 
to r^6 kpc (estimated from their placement of the inner 
4:1 circle in Figure 1, Plate 14, rather than with the dis- 
crepant value listed in Table 3 there). (Yet, we also note 
a four-armed appearance to the structure where spur-like 
features associated with either arm emerge, clear in CO 
emission and in the B-band). 



This result, namely that structure in HI shares the 
same pattern speed as the structure interior (exhibited 
solely by the molecular gas), admits a much stronger con- 
clusion than can be drawn from the CO, alone. Our rota- 
tion and resonance curves indicate that the inner pattern 
ends near its OLR (r^7 kpc; or possibly the outer 4:1 
resonance, within the errors) which would not be evident 
from the application of the TWR method to CO data 
covering r<3.5 kpc (corresponding to the central zone 
in Figures [1] and [5] with column density >10 20,6 cm -2 ). 
This finding is consistent with the theoretical expectation 
for the forbidden propagation of spiral density waves be- 
yond this resonance, w hich observations corroborate (see 
lElmegreen et al.lll989h . 

Likewise, corotation for this speed occurs at r=3.5 kpc, 
very near the transition between molecular and atomic 
gas (see Figure H]). With its lack of clear structure 
and relatively low surface density, this location in the 
gaseous disk seems consistent with an expected depopu- 
lation near CR owing to opposing torques (whereby gas 
is driven inward between ILR and CR and outward be- 
tween CR and OLR). (This impression, however, may be 
subject to the sensitivity of the CO map and the assumed 
X-i actor.) 

Interior to corotation, our determination of the disk 
angular rotation becomes less certain at a disadvan- 
tage to definitive resonance identifications. Neverthe- 
less, the two-armed spiral exhibited near r~2 kpc by 
the molecular gas-and which appears almost ring-like- 
is arguably located near the inner 4:1 resonance kpc, 
or UHR, where gas can accumulate on its pat h inward 
from corotation to ILR (|Buta fc Combes! Il996f ) . In ad- 
dition, as will be discussed in more detail in §3.2. 1[ the 
central concentration of gas in the form of the molecu- 
lar bar (with length ~1 kpc) appears to lie well within 
the CR. This is also true of the stell ar bar with a^2 
kpc (jKennev. Scoville fc Wilson! Il991h . which suggests 
a scenario quite different from observational findings 
in favor almost exclusively o f fast stellar bars (with 
l<r cfl /a;,<1.4; iCorsinil 120081 ). at least in early-type 
barred galaxies. Instead, the stellar bar here seems 
re miniscent of the slow b ars fo und i n the simulations 
of iCombes fc Elmegreenl ([1993) and iRautiainen et al.1 
((20081) . 

As for the pair of pattern speeds P 2 and P 3 out- 
side r t .i=6.7±0.25 kpc, we can similarly use the curves 
in Figure 0] to interpret their radial domains in terms of 
resonances. This is less straightforward, however, pri- 
marily because of our uncertainty in the poorly-fit ro- 
tation curve here, as well as the complex nature of the 
spiral structure visible in the HI surface density, which 
is characterized by po wer in several Fourier m odes (see 
Figure [5] and optically, Elm egreen et al.| [l9920. Notably, 
between r~6 kpc and 13 kpc, the structure visible in HI 
appears four-armed (see the top left panel in Figure [T] 
and Figure [5]). Meanwhile, near r=13 kpc a two-armed 
pattern starts to predominate again, as indicated in Fig- 
ure [T] and by Figure [5] where the power in the m=4 mode 
decreases relative to the m=2 mode. 

In this case, the transition r t .\ in Figure [4] would seem 
to associate the second speed fl p .2 with the four-armed 
spiral, here found to extend between the inner and outer 
4:1 resonances at r^6.8 kpc and r~14 kpc. (This is un- 
changed even with rotation curve models where the rise 



and fall in velocities between r=8 and 14 kpc in Fig- 
ure [3] are more closely (yet coarsely) fit, although the 
end of r2 p 2 more nearly approaches Q-k/2). The four- 
armed spiral in ESO 566-24 like wise extends between it s 
inner and outer 4:1 resonances (jRautiainen et alj|200l) . 
Our finding is also roughly consistent with the identifica- 
tion made by lElmegreen et "all |l992; in Figure 1, Panel 
14 there), namely that the zone dominated by the four- 
armed spiral is bounded by the inner 4:1 circle with r^6 
kpc and the CR circle with r~12 kpc. (We would argue, 
however, that the end of this zone occurs past the CR, 
which we find near r=10 kpc, depending on the rotation 
curve.) 

This is a scenario moreover favored by our TWR 
solution in that the transition to the four-armed spi- 
ral occurs at a resonance overlap; where P) i ends at 
OLR, within the uncertainties, f2 Pi 2 begins near its in- 
ner 4:1 resonance at r~7 kpc. Such an alignment of 
resonances, although not conclusively associated with 
mode-couplin g (i.e. as opposed to CR-ILR overlaps 
identified by iMasset fc Taggerlll997bD, has been id enti- 
fied in the simul a tions of IRautiainen fc Salol (|1999t ) and 
iDebattista et al.l (|2006h . as well as in the grand-design 
spiral M51 through application of the TWR method (Pa- 
per II). Presumably, this overlap is characteristic of a 
physical mechanism by which spiral s tructure can be sus- 
taine d over a large span in radius (Ra utiainen fc Salol 
fl999l) . 

The transition Tt,i between fl Pt 2 and Q p .3 may also be 
accompanied by resonance overlap: the outermost speed, 
which spans a radial domain well-matched to that of the 
two-armed spiral within 12<r<20 kpc (marked by the in- 
crease in to=2 power in Figured]), could begin at cither 
the ILR or inner 4:1 resonance. Large fractional errors 
on f2 Pi 3, of course, make this identification particularly 
vague. Likewise, the end of this pattern, as designated by 
the cut radius, can be only indefinitely related to realistic 
resonances (e.g. CR; at least with our current determi- 
nation of the angular rotation), although it is just be- 
yond the end of the ou ter two-armed spiral identified by 
lElmegreen et a l. (1992) near r=19 kpc (again, according 
to their Figure 1, Panel 14, rather than Table 3). In ad- 
dition, whether or not the resemblance to £1-k/2, H-k/4 
or fl suggested in Figure [4] is significant remains unclear 
at this point, especially in light of the nearly linear rise 
(and the asymmetry) in the rotation velocities outside 
r~19 kpc (not well modeled here; see Figure [3|) suggest- 
ing that the angular rotation curve may flatten out near 
this radius. The values of Q,-k/2 and JI-k/4 would be 
lower (and H+k/2 and fl+n/A higher) than in FigurelH 
in which case il p .3 may coincide with Q-k/A. 

Despite our uncertainty in the outermost speed, over- 
all, the TWR solution presents compelling evidence for 
extensive spiral structure described by multiple pattern 
speeds that are moreover related by their overlap at 
resonance. Given that this galaxy is tidally interact- 
ing (e.g. as argued bv lldetal l2002f ) . our pattern speed 
solution may be suggestiv e of the scenario explored by 
iSalo fc Laurika incn (20Q3) in simulations of M51. Those 
authors find that waves of higher and higher pattern 
speeds are excited as tidally-induced waves near Q-k/2 
propagate inward. But considering that resonance over- 
lap is not a feature in those simulations, our finding of a 
link between the speeds in M101 (and in M51; Paper II)- 
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similar to the overlaps so far de monstrated between bars 
and spirals in simulatio n (e.g. iMasset fc Taggerlll997bl : 
iRautiainen fc Saldll999l) -may imply that at least the in- 
ner disk is dominated by internal, rather than external, 
effects. The speeds and structure within r<r c ~22 kpc, 
however, have yet to be related to the asymmetric re- 
mainder of the disk. If the outer warp is a lop-sided, 
m=l mode, its pattern speed (not measured here) may 
be determined by, and reveal clues to, tidal encounter 
with the companion galaxies of M101. 

3.2.1. The inner pattern speed in M101 

As stated previously, the inner disk appears to sustain 
only a single constant pattern speed, Sip i=47±10 
km s _1 kpc _1 . This may be surprising, given the distinct 
bar pattern visible in the CO emission, which we might 
expect to end near CR and rotate with a pattern speed 
distinct from that of the spiral pattern immediately 
exterior. In order to determine whether our TWR 
solution is an accurate representation of the true speed, 
in this section we discuss the evidence in favor of only a 
single pattern speed in this zone, consider whether the 
measurement is an artifact of variation in the A-factor, 
and make comparisons to findings in the literature. 

For this investigation we consider the CO data alone. 
As might be expected given the radial domains of the 
two pattern speeds, we find that solutions calculated 
from the CO and HI maps, individually, supply pattern 
speed estimates nearly identical to those in either the 
inner or outer zones of the solution in FigurelH the TWR 
result from the CO maps, Sl p =48±8 km s _1 kpc^ 1 , is 
within the errors of the inner speed measured from the 
combined maps. Furthermore, the contribution made 
by the bar and inner spiral pattern in a given slice can 
be more easily identified when slices do not also contain 
information from throughout zone of the HI. (Note 
that, without the need for a cut radius, TWR solutions 
discussed below are all completely regularized.) 

Of course, the notably low sensitivity of the CO map- 
with the inter-arm regions and the zone surrounding the 
bar showing little , if any, emission (Figure [5J see also 
iHelfer et all l2003T) -itself raises significant doubts about 
whether the pattern speed can be accurately measured. 
According to Paper I, however, as it is engaged with the 
TWR calculation, regularization essentially attributes 
information within a given element of the TWR kernel 
to all elements at the same radius in all other slices, 
thereby compensating for zones which lack clear sig- 
natures of patterns. (Whether or not this is accurate 
depends primarily on how well the PA and inclination 
have been determined.) Errors that may be expected to 
arise in such low S/N regions are furthermore smoothly 
redistributed throughout the solution. 

Given that the net result is an effective extraction of 
the available information, we argue that our regularized 
solution is a reliable indication that the data do not 
support the measurement of more than a single pattern 
speed. In fact, when we impose a transition near the 
expected molecular bar end, solutions with 0.75<rt<1.6 
kpc either indicate no difference between the inner and 
outer pattern speeds, or the inner speed is only slightly 
higher than, and within the errors of, the outer. 

As evidence of the incompatibility of a distinct, higher 
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Fig. 7. — Plot of solution-reproduced (dots and crosses) and 
actual (open circles) integrals <v>i=b{/ f Hdx (see Paper II) as 
a function of slice position y at PA=42° calculated from BIMA 
SONG CO data for M101. The values associated with the best-fit 
single pattern speed solution (black dots) are plotted along with 
those of the two mock pattern speeds, f2 p ,i=115 km s — 1 kpc — 1 out 
to rt=X.l kpc and Q P; 2=48 km s _1 kpc _1 throughout the remain- 
der of CO-traced disk (red crosses). 

pattern speed with the data, in Figure [7] we show a 
comparison between the <u> (defined in Paper II and 
references therein) reproduced by the best-fit, constant 
pattern speed solution and by two mock pattern speeds: 
one for the bar, Sl&=115 km s _1 kpc _1 (chosen in order 
that the bar ends near CR) and one for the spiral, 
Sl s =48 km s~ kpc - (adopted from the value of the 
best-fit solution). There, the best-fit solution yields a 
significantly better fit to the data than the two speeds 
for slices inside |y|~1.0 kpc. 

This conclusion does not appear to depend on the 
PA or inclination, unlike in M51 where the radial varia- 
tion in the best-fit solution varies from PA to PA (Paper 
II). For all angle combinations a single, constant pattern 
speed yields a significantly better fit to the data than 
two, distinct speeds. Nor is the measurement sensitive 
to spatial variation in the CO-H2 conversion factor. 
Based on the metallicity gradient - 0.028±0.01 dex kpc 1 
measured from oxygen abundances (|Cedres et alj [2004). 
we modeled an increase in X with radius across the CO 
e mitting disk (ac cording to the galaxy-to-galaxy scaling 
of lArimoto et al.l fT996). In this case a single, constant 
pattern speed is not only once again the best-fit to the 
data, but at Sl p =53±6 km s~ kpc -1 , there is virtually 
no change to the measured value. 

Our result, namely, that there is no distinct pattern 
speed over the length of the molecular bar, therefore 
seems to be authentic. Observations of HII regions 
down stream of the molecular spiral arms (| Waller et al.l 
119971 ) also support our finding that the spiral arms lie in- 
side corotation (although our spiral speed i s much higher 
than the value inferred by I Waller et al.l (|1997h from 
measurements of the difference in position of the H Q and 
CO emission). The relative spati al offset betw e en th e 
CO and UV emission in Figure 3 of I Waller et"ai1 (|1997j ). 
while vague in places, seems sustained throughout the 
zone of the spiral and diminishes roughly near our CR 
radius. Since we find that the bar and spiral have the 
same pattern speed, this further implies that the bar 
ends well within its corotation radius. 
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Fig. 8. — Best-fit regularized solution for IC 342 with 
r c =9.8±0.19 kpc for PA=42°±3°. For this solution, bins exte- 
rior to r c (not shown) have been calculated without regulariza- 
tion. Dashed red lines and green horizontal error bars represent 
the dispersion in the pattern speeds and in rt and r c in solu- 
tions derived with a two-pattern speed model at the three PAs. 
The values in the zone of the bright spiral structure correspond to 
C p ,i=38±7 km s _1 kpc _1 out to r t =5.7±0.71 kpc and fi Pl 2=ll±6 
km s — 1 kpc — 1 out to r c . Curves for Q, f2±re/2 and C±k/4 are 
shown in black, cyan and blue. 

This scenario has bee n suggested by 
iKennev. Scoville fe Wilson! (|1991h who, in first re- 
porting on the 25° offset between the molecular and 
stellar bar position angles, no ted a resemblance to th e 
hydrodynamical simulations of iCombes fc Ger in (1985). 
Such a decoupled central gas concentration is found 
when the stellar bar is slow and drives spirals that 
develop inside corotation, as we suggest here. 

A slow bar, ending well inside corotation, also seems to 
be reconciled with central DM content in M101 implied 
in the multi-component simulations of Petitp as et al.l 
( 2004), which best reproduce the observed velocity field 
with a minimum disk. Where the DM contribution to the 
gravit ational potential is large, iDebattista fc Sellwoodl 
(2000) find that interaction between a bar and the DM 
halo through dy n amical friction decelerates t he bar 
(jWeinberd 1198a IDebattista fc Sellwoodl |2000D . such 
that it grows in length disproportional to the greater 
increase in Tcr- It should be noted, however, that the 
slowness of the bar implied in this case is a matt er of 
debate (cf. ISellwoodl 120081 iDubinski et "all (|2009t) and 
I Weinberg fc Katd 120071 who argue that slow bars are an 
artifact of the resolution of the N-body simulations used 
to model the interaction). 

3.3. IC 342 

The TWR calculation for this galaxy proceeds with the 
use of a radial bin width Ar=0.41 kpc (as established by 
the resolution of the combined maps) using slices out to 
1 2/ = 1 3 . 3 kpc. The best-fit solution, plotted in Figure 
[S] together with angular rotation and resonance curves, 
measures two distinct pattern speeds inside the cut ra- 
dius r c =9.8±0.19 kpc. As in M101, we find that this cut 
radius accurately reflects the location where the distor- 
tion/warp in the outer disk begins, e. g. as identified in 
the H I surface density/velocity field (jCrosthwaite et all 
120001 see Figure [TJ. 

Furthermore, the speed r2 p i=38±7 km s _1 kpc _1 in- 
side r t =5.7±0.7 kpc is nearly identical to the value 
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Fig. 9.— Total H 2 +HI surface density in IC 342 highlighting the 
structure inside ~9.8 kpc. The transition rt marking the extent 
of the inner pattern speed in our solution is shown as a dashed 
black circle. The horizontal bar near the bottom left indicates the 
physical scale. 

iCrosthwaite et all (|2001h estimate based on their deter- 
mination rcij=4.2±0.7 kpc from inspection of the field of 
velocity residuals. According to the solution in Figure 
as parameterized by r t this pattern terminates at either 
corotation (here at r^5 kpc; conditio nal on the rotation 
curve , which is different from that of ICrosthwaite et al.l 
2001) or OLR (within the uncertainties) in favor of the 
second, lower speed ri Pi2 =ll±6 km s _1 kpc _1 . 

The transition between the two speeds occurs very near 
where the spiral structure becomes visibly four-armed 
(see Figure H]) . This seems to establish that the outer 
arms are, in fact, best described by a speed that is dis- 
tinct fro m that of the spiral interi or, as previously sug- 
gested by ICrosthwaite et all (|2000f ). Moreover, Figure [5] 
suggests a possible resonant link between the two speeds 
in IC 342, in that the lower speed begins near its ILR 
(or possibly its inner 4:1 resonance). The low speed 
also seems to end at the 4:1 resonance, and possibly CR 
(within the uncertainties). However, the cut radius r c 
in Figure [5] seems less an accurate bound on the four- 
armed spiral, which continues out to r~ll kpc (see the 
zeroth- moment map in Figure [T]), than an indication of 
where the warp begins. In this case, whether the warp 
and four-armed spiral are related or share a common ori- 
gin is unclear. 

Although we might also infer from Figure [8] that the 
inner speed lacks an ILR, our uncertainty in the rota- 
tion curve at the center is high, given the low resolu- 
tion of the maps. Without more confidence in the be- 
havior of the £1-k/2 curve, in particular, it remains un- 
clear whether this speed, and the structure it describes, 
has an inner bound. Currently, then, both the spiral 
structure tra ced by CO and the distin ct bar with length 
as ~1.5 kpc (jCrosthwaite et all l2000h appear to rotate 
with speed f7 p ~40 km s _1 kpc _1 . As in M101, this no- 
tably implies that the bar ends well inside its corotation 
radius, a circumstance which may also be evidenced by 
the ^9° offset between the molec ular and stellar bar ma- 
jor axes (jCrosthwaite et af1l2001h . as similarly discussed 
in 33~2~i1 
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Fig. 10. — Best-fit regularized solutions for NGC 3938 with 
PA=21°±2°. Dashed red lines for both solutions, fl£ (thick line) 

and O? (thin line) represent the dispersion in the pattern speeds 
at the three PAs. The green horizontal error bar shows the disper- 
sion from PA to PA in rt in solutions derived with a two-pattern 
speed model. In the two-speed solution, the value in the inner zone 
corresponds to £7^=53^19.2 km s _1 kpc _1 out to rt=3.4±0.7 kpc 

followed by £2^ 2 = 7.7±1.6 km s — 1 kpc — 1 , while the single, constant 
speed solution measures n^=6.6±2.1 km s _1 kpc _1 . Curves for fl 
and Q±k/2 are shown in black and cyan. 

On the other hand, the molecular bar and inner spi- 
ral could have different pattern speeds measurable, in 
principle, but for the width of the radial bins; informa- 
tion from the zone of the molecular gas, which covers a 
relatively small extent, is limited to only a minor con- 
tribution in solutions. Although the two-pattern speed 
solution in Figure [5] yields a significantly better fit to 
the data, a higher pattern speed can, in fact, be recog- 
nized in solutions modeled with three distinct pattern 
speeds: the lowest xt solution of this type at PA=42° 
favors a transition at rj=1.6±0.45 kpc from f2 p =34±9 
km s _1 kpc _1 to a higher fl l p —70±12 km s _1 kpc _1 , with 

Sl p =12±4 km s _1 kpc _1 between r t =5.3±0.64 kpc and 
r c =9.8±0.19 kpc. (Errors represent the dispersion in 
three-speed solutions where the inner-most transition, 
found best at the optimal PA, is held fixed from PA to 
PA). This speed is particularly compelling since the bar 
would end very near CR, and at this radius overlap with 
the lower, spiral pattern's inner 4:1 resonance. 

However, we emphasize that this identification is in- 
conclusive at this point; the zone of the higher pattern 
speed is covered by only four radial bins, or less than 
11% of the disk. (Note that the CO-which has a much 
higher resolution than the Hl-alone cannot be used as 
a continuity-obeying kinematic tracer, unlike in M101.) 
Sensitivity to variation in the CO-to-Eb conversion factor 
(e.g. giv en the metallicity gradient in this galaxy mea- 
sured bv lVila-Costas fc Edm unds 1993) is also currently 
untestable; with such a large radial bin width any varia- 
tion in X that we might accommodate would be crude at 
best and, with only four bins covering the CO emission, 
would likely yield an insignificant result. 

3.4. NGC 3938 

Given the resolution of the HI data (to which the 
BIMA CO cube had been smoothed prior to the com- 
bination) and, as a consequence, the rather large bin 
width Ar=1.7 kpc, our TWR solutions calculated in 
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Fig. 11.— Total H 2 +HI surface density in NGC 3938 highlighting 
the structure inside ~15 kpc. The transition rt marking the extent 
of the inner pattern speed in solution Cl^ is shown as a dashed 
black circle. The horizontal bar near the bottom left indicates the 
physical scale. 

slices | 1 < 13. 1 kpc contain only 8 radial bins. This, to- 
gether with the low inclination of the disk-which com- 
promises the signatures of departures from axisymmetry 
in the velocity field-limits how well radial variation in 
the pattern speed can be detected, if present. Under 
this scenario, we find that two TWR solutions calcu- 
lated with regularization over all radial bins (i.e. with 
no r c ) fit the data equally well. The first ascribes a sin- 
gle, constant pattern speed Qp =6.6±2.1 km s _1 kpc _1 
to structure throughout the disk, while the second in- 
corporates the measurement of a distinct, higher pattern 
speed r^ 1 =53±9.2 km s _1 kpc _1 inside r t =3A kpc (in 

addition to the lower f^ 2 =7.7±1.6 km s _1 kpc _1 ). 

The measurement of any pattern speed at all may be 
surprising, given the low inclination angle and the ap- 
parent lack of strong, organized structure in both the HI 
and (smoothed) CO emission analyzed here. Compar- 
ison with the rotation curves in Figure [10] greatly aids 
our interpretation of the measurements, although any 
conclusions that might be drawn are tenuous at best; for 
this low inclination, rotation curves are highly uncertain 
(though similar to those thr ough out the literature, e.g. 
[ Jimenez- Vicente et al.l 119991 and ICombes fc Becauaertl 
1997). 

With this in mind, we find it noteworthy that the low 
speed fiy (or $lp 2 ) IS vei T near in the outer disk 

(Figure [ini an d possibly near O-k/4, as well). This sim- 
ilarity seems reminiscent of a pattern formed by closed, 
precessing elliptical orbits guided by the near-constancy 
of Q-k/2. But since there is no dominant m=2 Fourier 
component visible in the 21 cm emission or in th e NIR (as 
observed by ICastro- Rodriguez fc Garzonl 120031 ) at these 
radii, it is unclear that this speed is associated with spi- 
ral structure, at all. A lack of spiral streaming motions 
in the atomic gas ( as also observed in the ioniz ed com- 
ponent of the ISM; Uimenez- Vicente et al.|[T999l ). would 
furthermore lead us to expect that, for any asymmetry 
in the surface density, a value much closer to the disk 
angular rotation frequency f2 would be measured. 



12 



At radii r<6 kpc, on the other hand, clear spiral struc- 
tu re is identifiable in the NIR (see the K-band image 
in ICastro- Rodriguez fc Garzon| [2003). At the very least, 
this seems to suggest that the two-pattern speed solution 
flp (with r t as marked in Figure fTT|) may be more realis- 
tic than f2p ; the faster inner speed seems more compati- 
ble with structure at these radii than O p ~7 km s _1 kpc _1 
since, even with uncertainty in f2 at the inner radii and 
the near-constancy of CI-k/2 at the outer radii, the low 
pattern speed's ILR would lie at r>5 kpc. 

Nonetheless, whether f2^j is a reliable measurement of 
the expected higher pattern speed is difficult to assert. 
The transition r t (see Figure fTT)) lies just inside the edge 
of the CO emission tracing a two-ar med spiral (clear i n 
the unsmoothed BIMA SONG map; iHelfer et al.l [2003h . 
and so the measurement presumably reflects the speed of 
this pattern. However, even though this speed appears 
to cover a quarter of the disk, in general we expect that 
two radial bins would not supply sufficient leverage on 
the measurement. (This is arguably the reason why the 
two TWR solutions are indistinguishable). Indeed, with 
so few bins in this zone we are prevented from diagnosing 
radial variation of order higher than zero. Furthermore, 
the accuracy of the measurement, if real, depends on how 
accurately information about the two speeds can be sepa- 
rated according to the parameterization of the transition 
r t —3A kpc, which here may be severely limited by the 
size of the radial bin width. (Note that, as in IC 342, the 
unsmoothed 6" resolution CO alone cannot be used as 
a continuity-obeying kinematic tracer, since ISM is not 
molecule-dominated over the extent of the CO emission. 
Tests for sensitivity to variation in the CO-to-H2 conver- 
sion factor are also not well-accommodated with so few, 
large radial bins). 

These shortcomings notwithstanding, the distinct mea- 
surement inside r~3 kpc seems credible if only that 
it is consistent with values indicated by two indepen- 
dent appro aches using diffe r ent t racers of the spiral 
structure. iKorchagin et al.l (|2005f ) find with a global 
modal approach that the dominant m=3 and m=4 
modes in their multi-wavelength observation-based simu- 
lations are well-described with a pattern speed near ^55 
km s _1 kpc _1 . iMartinez-Garcia et al . (2009) propose a 
similar pattern speed, ^47 k m s _1 kpc _1 (adopted from 
IMartinez-Garcia et al]|2009l for D=11.3 Mpc), based on 
the azimuthal age gradient across the spiral arms calcu- 
lated from optical images. 

From Figure [TU] we might also argue that, within the 
uncertainties, ends near corotation at rcR—2.b kpc 
and so represents a physically realistic scenario, also re- 
cently identified in M51 (Paper II). However, this speed 
might just as reasonably end at OLR, depending on the 
rotation curve, the determination of rt , and the measure- 
ment itself. Also, while this pattern appears to lack an 
ILR, this is uncertain given that the angular velocities are 
the least well-determined at the inner radii. At this point 
there also does not seem to be a clear relation between 
the two pattern speeds (e.g. overlapping CR and ILR; 
iMasset fe Taggerlll997bD . especially since resonances for 
the outer speed f2^ 2 are difficult to establish, given the 
flatness of the resonance curves. 

While our confidence in the TWR measurement is tem- 



pered by the low inclination of the disk, the value of this 
analysis lies in the prospect of renewed perspective on 
the nature of the structure in this galaxy. Although in 
no way definitive (based on the TWR method, alone), 
the lower measurement near 51-k/2 (either or f^ 2 )> 
in particular, may be compatible with several pieces of 
evidence that po int to the influence of the DM h alo first 
contemplated by Ivan der Kruit fc Shostak (1982). 

As investigated bv lFrenk et al.l (|1988f) or bv lJod (|1997H . 
for example, the DM halo may be responsible for the ap- 
pearance of structure in this isolated galaxy, which other- 
wise seems difficult to reconcile with the finding that the 
gas disk is everywhere subcr itical to gravitational insta- 
bility (as remarked upon bv lCombes fe Becauaertlll997l . 
and aside from an alternative increase in instability pos- 
sible through the c oupling be tween multiple components 
in the system, e.g. Uog||1992T ). 

Specifically, w here the Toom r e stab ility parameter in 
the g as is high, iBureau et alJ (|1999t ) and iFrenk et al.l 
(1988) suggest that the torque due to a triaxial halo 
(predicted by Cold Dark Matter simulations of hierar- 
chica l structure formation; e.g. iDubinski fc Carlbergj 
fl99l[ ) with slow fi gure rotation can drive structure in 
extended HI disks. Material (though not exclusive) to 
this scenario, the disk of NGC 3938 exhibits intrin- 
sic ellipticity: as sim ilarly diagnosed in ot her eccen- 
tric nearby spirals by lAndersen et alJ ([2001), the pho- 
tometric and kinematic position angles of NGC 3938 
are offset by nearly 50° (iDaigle ef aj] |2006|) . In fact, 
iCastro-Rodriguez fe Garzonl 120031 measure an ellipticity 
e,/=0.11 in the J-band. So while the ellipticity may have 
arisen with an asymmetric accretion of matter, for ex- 
ample, it could also reflect an m—2 potential pertur- 
bati on occasioned by tria xiality in the DM halo (see, 
e.g., lAndersen et al.l [20011 and references therein). In 
this case, the TWR solution may represent a measure of 
the rotation speed of such a halo, if the structure in the 
out er disk (which, extends further i n HI than in the opti- 
cal; Ivan der Kruit fc Shostaklll 982) arises in the manner 
considered bv lFrenk et al.l (|1988l ). Currently, however, it 
is not clear that the structure in NGC 3938 is compatible 
with th e strong two-armed patterns found in the simula- 
tions of lFrenk et all l| 1988ft . 

Alternatively, the overall instability of the disk might 
be increased by a global mass asymmetry resulting from 
th e rn=l perturbat ion t o the halo potential considered 
by [Joi (|1997l 120001 and l2002h . which is also a possible 
source of the ellipticity. This latter scenario seems to 
be compatible with other observa ble c haracteristics of 
this g alaxy: like iTullv et al.l (|1996h and iBournaud et all 
(2005), we find the HI disk to be slightly lop-sided toward 
the north, and this m—1 asymmetry appears in the HI 
velocity field, as well. ( As previously dem onstrated with 
the Ha observations of lDaigle et al.ll2006l . we measure a 
^10-20 km s _1 difference between the approaching and 
receding sides beyond r^ll kpc.) In this case, the TWR 
measurement may reflect the response of the gas (and 
stars) to this pertu rbation if the elliptical orbits calcu- 
lated by Pod (|1997ft precess with frequency VL-k/2. 

Again, our TWR measurement does not confirm, or 
distinguish between, these possibilities. In either of the 
two scenarios it is also not clear if, or how, the inner, 
higher speed may relate to the lower speed. Higher reso- 
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Fig. 12.— Best-fit regularized solution for NGC 3344 with 
r c =6.8±0.43 kpc for PA=155°±2°. For this solution, bins exte- 
rior to r c (not shown) have been calculated without regularization. 
Dashed red lines and the green horizontal error bar represent the 
dispersion in the pattern speeds and r c in solutions derived at the 
three PAs. The value in the zone of the bright spiral structure 
corresponds to f2 Pi i=44±4 km s kpc out to r c . Curves for O, 
f2±ft/2 and £7±k/4 are shown in black, cyan and blue. 

lution observations are necessary to first establish which 
of the two solutions, Vl^ or £1^ , is the most appropriate 
for this galaxy. Only then will it be possible to more 
rigorously explore the relation between the two pattern 
speeds suggested by the solution £1^ and perhaps then 
recognize the true nature of the low speed. 

3.5. NGC 3344 

The best-fit solution calculated in slices out to 
|?/| = 19.4 kpc with radial bin width Ar=0.75 kpc (cor- 
responding to the ~20" resolution of the HI cube) is 
shown in Figure [12] along with rotation and resonance 
curves. Immediately we notice that, according to our 
determination of the cut radius r c =6.8±0.43 kpc-near 
the location where the warp has been identi fied to begin 
(r~7 kpc; IVerdes-Montenegro et al.l l2000f )— this spiral 
pattern ends at OLR, within the uncertainties. The 
inner bound for the single, constant speed 17 p =44±4 

km s _1 kpc _1 in Figure [12] on the other hand, is less 
clearly well established here. The paucity of 21 cm 
emission inside r~l kpc not only leads us to suspect 
that the measurement may not be valid all the way to 
the center, as suggested by our solution in Figure [12] 
(calculated with bins covering all radii r>0), but it also 
prevents us from constraining the angular rotation curve 
there (and the position of the ILR) with confidence. 

Even so, according to our estimate r/i/j^l.O kpc, the 
spiral could reasonably begin at the inner ring, as iden - 
tified in the optical bv IVerdes-Montenegro et al.l (2000). 
In this case, with the single pattern speed implied by 
our TWR solution, bo th the inner and outer rings, a t 
r=l kpc and r=7 kpc (jVerdes- Montenegro et al.l 12000). 
are arguably associated with spiral's inner and outer 
Lindblad resonances. We note that this depends in no 
small part on the TWR method; the separation achieved 
by the cut radius between the strong spiral structure 
and the m=l distortion grants a measurement for the 
former constant speed which is significantly different 
from that implied by the traditional TW method, 
fi p =26±6 km s _1 kpc _1 (calculated with the same slices, 
kinematic parameters, and limits of integration as in the 



TWR calculation). 

The model-independent TWR method, more- 
over, resolves the ambiguity in the study by 
IVerdes-Montenegro et al.l ([2000). There, the uncer- 
tainty in the rotation curve precluded the reliable 
association of the inner and outer rings with res- 
onances, and the estimation of a pattern speed 
(|Verdes-Montenegro et~aT1 12000D . With the identifica- 
tion above we can confirm that the outer pseudoring at 
r^7 kpc o riginates with the spiral, as op posed to the 
bar, which IVerdes-Montenegro et al.l (|2000f ) argue is too 
small to have an OLR at this radius. 

The inner ring, on the other hand, seem more 
likely related to the bar (which it surrounds; 
IVerdes-Montenegro et all [2000), given the proxim- 
ity of these two features. That we find the spiral's ILR 
nearly coincident with the location of the inner ring 
might then be a manifestation of resonance overlap 
between the speeds of the bar and spiral. If the spiral 
with lower speed is driven at resonance, for example, our 
TWR solution may suggest a reasonable speed for the 
bar, which we are prevented from measuring here. In 
particular, if we associate the inner, ILR ring with inner 
extent of the spiral (inside of which the HI density falls 
off), such a scenario could be realistically achieved via 
CR-ILR mode-coupling, where the bar rotates with pat- 
tern speed Q P >150 km s~ 1 kpc ~ 1 for rqR^l.Q kpc. The 
bar, with length a_e~0.7 kpc (|Verdes-Montenegro et al.l 
2000), in this case would be reasonably fast, with 
rcij/as^l-4. Alternatively, we might expect a bar 
pattern speed il p >80 km s~ 1 kpc _1 , in the event of the 
CR-inner 4:1 resonance overlap discussed in £13.21 in 
which case rc_R/as~3.5. In the former (latter) scenario, 
the inner ring could be located near the bar's CR (inner 
4:1 resonance), at least with the rotation and resonance 
curves derived here. 

Our TWR solution also se ems to endorse the sce- 
nario speculated upon by IVerdes-Montenegro et al.l 
(2000) for the origin of the warp in this galaxy, 
namely through t he no n-linear coupling proposed by 
iMasset &; Tagged (]1997al) . In this scheme, the warp 
(manifest in the lop-sided part of the HI distribution, 
which extends 34% further to ward the SE than NW; 
IVerdes-Montenegro et alJfeOOOh would originate through 
a coupling at OLR between a spiral density wave and 
two "warp waves" . Since we have established that the 
spiral structure ends at OLR, in the ab sence of tidal 
inter action with a companion galaxy (|Thilker et al.l 
l2007f ). the spiral pattern itself could be capable of 
generating the warp in the outer disk. 

4. SUMMARY AND CONCLUSIONS 

Direct pattern speed measurement affords an observa- 
tional resolution for several fundamental issues in the na- 
ture and origin of spirals. The relation between bar and 
spiral pattern speeds, the number and radial domains of 
patterns speeds that can be sustained in the disk, and 
whether spiral structure is steady or winding, for ex- 
ample, can all be established with knowledge of spiral 
speeds and their radial variation. The TWR method, a 
technique for measuring radially varying pattern speeds, 
supplies us with the first such measurements to address 
these issues. 

In this paper we have applied the TWR method to ob- 
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servations of CO and HI in four spiral galaxies. For this 
work, we have expanded the number of spiral galaxies to 
which TW-type calculations can be applied by consider- 
ing, for the first time, the combination of both molecular 
and atomic gas. Together, the two ISM phases better 
meet the continuity requirement of the method and also 
afford greater insight into whether (and how) multiple 
spiral pattern speeds extend over a large range of radii. 
This has notably increased the sample size of galaxies 
analyzed with the TWR method so far. 

Our TWR solutions for M101 (S3), which, of all the 
solutions presented here, are equipped with the small- 
est radial bins, very clearly show radial variation in the 
pattern speed across the disk. Within the inner 3', we 
find convincing evidence that the bar and spiral have the 
same pattern speed (distinct from an exterior speed), and 
that both lie inside corotation. This represents a scenario 
quite different from recent findings in fav or almost ex- 
clusively of fast bars (with 1 <m r I ah < 1 . 4; ICorsinill2008f) 
albeit in mostly early-type barred galaxies. 

In MIOl-as well as in IC 342 and possibly NGC 3938- 
we also find that the extensive spiral structure there is 
best described with more than a single pattern speed. 
Furthermore, in both M101 and IC 342 we find evidence 
that the transition between two speeds coincides with 
resonance overlap. In the former, the transition between 
the two- and four-armed patterns occurs near the overlap 
of the inner's OLR and the outer's inner 4:1 resonance. 
In the latter case, the inner pattern transitions to a four- 
armed spiral at an OLR-ILR overlap. 

Although the specific resonance overlaps identified 
in M101 and IC 342 have not been conclusively 
demonstrated as t r ue instances of mode-cou pling (e.g. 
ISvgnet et all 119881 . iMasset fc Tagged Il997bf) . together 
with the CR-inner 4:1 resonance overlap in M51 (Pa- 
per II), this work suggests that there exist several pos- 
sibilities dictated by resonance overlap by which exten- 
sive spiral structure can be sustained. Our findings are 
qualitatively similar to t he barred spiral simulations of 
iRautiainen fc S alo (1999|), but we note that the spatially 
coincident existence of multiple modes with different pat- 
tern speeds, as often found there beyond the bar, is 
untested with the TWR calculation here. Even so, the 
outer four-armed spirals in both M101 and IC 342 are 
arguably excited at resonance by the patterns interior. 

On the other hand, in our sample we also find exten- 
sive spiral structure in the absence of clear resonance 
overlap. The speed characteristic of the outer, flocculent 
structure in NGC 3938, for instance, does not seem to be 
clearly related to that of the inner, two-armed spiral, if 
the two speeds are distinct. In addition, the tight spiral 
structure throughout the HI emitting disk in NGC 3344 
is best described with a single, constant pattern speed. 

In no case do we find that pattern speed is a smoothly 
decreasing function of radius, as might be expected for a 
winding spiral. Our cut procedure, however, may intro- 
duce a bias against such spirals if they exist preferentially 
in the outer regions of the disks in our sample; the ra- 
dial bins covering the outer portion of three of the four 
galaxies where a warp is evident are excluded from our 



models and calculated without regularization. 

Even though in this case we cannot characterize the 
patterns that may distinguish the warped regions of the 
disks in our sample, accurate measurements for the pat- 
tern speeds of the structure interior can, themselves, rec- 
ommend different mechanisms by which warps are ex- 
pected to originate and evolve. The galaxies in this 
sample are consistent with various existing explanations. 
Where the spiral structure ends at OLR in a warped, 
isolated galaxy (e.g. NGC 3344), for example, spirals 
themselves may be capable of exciting the warp through 
the no n-linear coupling proposed by IMasset fc Taggerl 
(|1997a| L In the absence of both a companion and this 
resonance boundary, but where the disk is intrinsically el- 
liptical, as in NGC 3938, the ellipticity, the warp and/or 
the instabilities in the disk may arise with a DM halo or 
the asymmetric accretion of matter (discussed in § 13. 4[) . 
An obvious companion, on the other hand, is a clear, 
potential source of tidal perturbation to the disk (e.g. 
M101 and IC 342). Although an accompanying warp or 
distortion may not necessarily relate to the spiral struc- 
ture at resonance, where the speed of outer structure is 
near the Q-n/2 curve, for example, this structure may 
be reson antly excited (at ILR) by the external perturber 
(see e.g. ISalo fc Laurika incn 20Q3). 

These issues, and the existence of extensive spiral 
structure, are best developed with a larger sample of 
galaxies, where meaningful trends will be more clearly 
established. For greatest effect, the TWR calculation 
should be applied to higher resolution HI observations 
than we analyze here. In all of the galaxies in our sam- 
ple (excepting M101) the resolution of the HI data limits 
our confidence in how well radial variation is constrained. 
It also inhibits our leverage on information extracted in 
the zone traced by CO emission (which we smooth to the 
resolution of the HI in combining the two observations). 

A large radial bin width likewise prevents us from sat- 
isfactorily testing the effect of variation in the CO-to- 
H2 conversion factor on TWR pattern speed estimates. 
Future studies with higher resolution observations (and 
therefore more radial bins across the CO-emitting disk) 
may be able to better test for the effects of a gradient 
or arm-interarm variations in X. We note, though, that 
these types of variations we re found to have a negligi- 
ble e ffect on TW estimates (jZimmer. Rand fc McGrawl 
2004), and here we find that low- level variation in the 
A"-factor produces very little change in our TWR solu- 
tions for M101. 

In general, these types of TWR measurements in 
nearby galaxies are invaluable for interpreting observa- 
tional studies of the evolution of bar and disk parameters 
now possible with HST. They also promise to be signifi- 
cant for studies at intermediate redshift, and future stud- 
ies with JWST, which will extend structural parameter 
measurement to larger redshifts and allow smaller bars 
and the earlier evolution of disks to be studied. 

We would like to thank Pertti Rautiainen for his con- 
structive comments on this manuscript. This material is 
based on work partially supported by the National Sci- 
ence Foundation under grant AST 03-06958 to R. J. R. 
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